The Current-Temperature Phase Diagram of Layered Superconductors 
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The behavior of clean layered superconductors in the presence of a finite electric current and in zero- 
magnetic field behavior is addressed. The structure of the current temperature phase diagram and 
the properties of each of the four regions will be explained. We will discuss the expected current 
voltage and resistance characteristics of each region as well as the effects of finite size and weak 
disorder on the phase diagram. In addition, the reason for which a weakly non-ohmic region exists 
above the transition temperature will be explained. 
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I. INTRODUCTION 

The understanding of a current's effect on the behavior 
of vortices in clean layered superconductors is important 
on it owns merit and is essential for the interpretation of 
measurements which use an electric current as a probe. 
This issue has been addressed by many authors. Glaz- 
man and Koshelev [0 found that in contrast to quasi-two- 
dimensional superconducting films, there is a non-zero 
critical current because a minimum current is needed to 
overcome the attraction due to the Josephson coupling 
between the layers and to drag apart vortex pairs. Jensen 
and Minnhagen 0] have generalized the two-dimensional 
current voltage (I-V) relation to the layered case by con- 
sidering the effect of the Josephson coupling on the in- 
tralayer vortex attraction. They find that the voltage 
goes from being exponentially small to finite with the 
following nonlinear dependence on current: 



V = k(T)I(I - ll{T)) 



a(T) 



(1) 



where ij(T') represents the temperature dependent 
threshold current needed to overcome the Josephson at- 
traction and k(T) and a(T) are generally taken to be 
independent of current. The relevance of this equation 
to such layered materials as the high-temperature super- 
conductors (HTSC's) has been established by several ex- 
perimental groups. [[|-|| 

A rigorous self-consistent study of the effects of a cur- 
rent on the critical behavior of vortices in a clean lay- 
ered system was conducted by the author using a renor- 
malization group analysis which culminated in the I-T 
phase diagram. H As will be described below, the I-T 
space was found to be divided into four regions separated 
by three characteristic currents which included a second 
order phase transition line. This too has been shown 
to describe well the behavior of the copper oxides. |-|] 
Martynovich and Artemov |l(J have also studied this sys- 
tem in zero and finite current and also find a second-order 
phase transition, but only at larger currents. Langevin 
simulations have also been used to study the effect of a 
driving current on an anisotropic Josephson junction ar- 
ray, [lilt] Gr0nbech-Jensen, Dominguez, and Bishop con- 



firmed the structure of the I-T phase diagram || finding 
that the phase transition occurs at a temperature higher 
than the temperature marking the onset of resistance and 
that the slopes of two of the characteristic currents (I^(T) 
and I C (T)) were linear. While the theoretical and simula- 
tional work seems sufficient for describing clean, layered 
systems, there is now evidence that finite size effects can 
produce a linear current dependence in the V-I relations 
at small current which give way to a nonlinear behavior 
at larger currents. fl|l^-|l4| This has also been seen in 
2D systems. [Jl5],[l6) These and other properties of lay- 
ered superconductors in the presence of a current need 
to be addressed. 

In this paper, we explain in more detail and expand 
upon the results of Ref. |^. In particular the details 
of how the I-T phase diagram was arrived at (Section 
[IB), the nature of the actual phase transition (Section 
UC), and the properties of each region (Section III) are 
explained more thoroughly. The origin of a weak non- 
ohmic contribution above the transition temperature will 



be discussed in Sections III C and [II D as will the effects 



of free vortices created due to finite size effects (Section 
[V). Finally, the effects of weak disorder as calculated 
via the replica technique will be reported (Section ^). 



II. MODEL AND DERIVATION OF I-T PHASE 
DIAGRAM 

A. Model and Recursion Relations 

We begin with a description of our model which has 
been described in detail elsewhere and so will be 

kept brief here. Vortices in layered systems are modeled 
as a "gas" of charges in a stack of weakly coupled lay- 
ers interacting via two-body interactions which approxi- 
mate those of vortices in layered superconductors. Inter- 
actions between vortices in the same layer and between 
vortices in neighboring layers are included but interac- 
tions between vortices separated by more than one layer 
are neglected. The effect of the current J is to exert a 
constant force on vortices in a direction perpendicular to 
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J but in opposite directions for oppositely charged vor- 
tices. Knowing the interaction potentials and the effect 
of the current, one can write down the partition function 
and perform a real-space renormalization group study on 
it in the fashion of Kosterlitz. [jl9| The results are, 

dx/de = 2y 2 [l — (1/16)A + J 2 ], (2) 
dy/de = 2y[x+ (l/2)Aln A]/(l + x) + Jy, (3) 
dA/de = 2A[1 - 4y 2 (I + J 2 )/(l + x)], (4) 

dJ/de = J. (5) 

A is the ratio of the interlayer coupling to the intralayer 
coupling p 2 , e — ln(7/£o) (where I is the relevant length 
scale), x = A/(f3p 2 ) — 1, and £ is the zero- temperature 
correlation length, y = exp(— f}E c )/l 2 is the fugacity, 
where E c is the "core energy." Note that the current 
affects not only the interaction strength parameters x 
and A but also the fugacity directly. This is because the 
current lowers the core energy of vortices and therefore 
directly affects the fugacity. 

B. The Correlation Length and the I-T phase 
diagram 

The I-T phase diagram is derived by studying the cor- 
relation length £(T) which is related to the value of I 
at which one terminates the integration of the recursion 
relations. Because the calculation leading to Eqs. (Q)- 
(|b|) assumed small vortex density (y), current, and inter- 
layer coupling, the integration of the recursion relations 
must be terminated when any of these quantities becomes 
large. The value of / at the cutoff is called l max and be- 
cause only one length (the correlation length) dominates 
the critical behavior, the following identification follows: 
lmax = The temperature dependence of £ is found by 
considering the first integral of the 2D recursion relations: 
y 2 + 2 ln(l + x) - 2x = c = y 2 + 2 ln(l + x l )-2x i , where 
the subscript i denotes the initial or "bare" value of that 
parameter, c is a constant of integration and it can be 
shown that c oc (T - T KT )/T KT @ where T KT is the 
2D transition temperature. Because we have assumed 
the interlayer coupling to be very weak in our model, it 
is expected that c should be linear in temperature in the 
immediate vicinity of Tkt for our system. 

Plotting ln(Z mQX /£o) versus c for various currents gives 
an indication of the temperature dependence of £(T) and 
the plots are qualitatively similar to that of Fig. 2 in 
Ref. pq . One finds that there is a peak in £(T) which 
occurs at the transition temperature. As one increases 
the current, the peak shifts to the left corresponding to a 



decrease in the transition temperature. This process al- 
lows one to determine T C (J) (or equivalently I C (T)) which 
decreases linearly with current. I C (T) is plotted in Fig. |l| 
as calculated for Xi = 0.5 and \ — 10 -6 . (j/j and Ji were 
varied to sweep through I-T space.) 
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FIG. 1. The I-T phase space with the three characteris- 
tic currents derived by studying the correlation lengths as ex- 
plained in the text. The x-axis label c is prop ortional to 
T/Tkt - 1. 

Two more characteristic currents can be derived by 
comparing the correlation length for the layered, finite- 
current case (referred to as the full correlation length) 
to the 2D finite-current correlation length and to the 
layered zero-current correlation length. In Fig. 1 of 
Ref. Q, the three correlation lengths are plotted. As 
expected, the 2D and full correlation lengths coincide 
at larger temperatures signifying the 2D behavior of the 
layered system for temperatures larger than T C (I). The 
small temperature range over which these two correlation 
lengths separate is taken to be the 3D/2D crossover re- 
gion and marks another characteristic temperature T 2 (J) 
(or equivalently I 2 (T)) which is also found to vary lin- 
early with current. (See Fig. |l|.) 

The third characteristic current comes from comparing 
the full correlation length with the layered, zero-current 
correlation length. Where these two correlation lengths 
begin to separate is the temperature at which the cur- 
rent starts to affect the system and identifies T}{I) (or 
ll{T)). Like I 2 {T) 1 it is a crossover current and not a 
phase transition. By varying the current, it is found that 
Ic(T) is linear at larger currents and is roughly parallel 
to I C (T) but crosses over to join I C (T) at small currents. 
In other words, ll{T c ) — I C (T C ) as illustrated in Fig. [j]. 
As we shall discuss below, what happens in finite current 
at two different temperatures, occurs at one temperature 
in zero current. 

While I C (T) can be determined to some precision, the 
other two characteristic currents cannot be since they 
are crossover currents. One can determine them in dif- 
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ferent ways which yield results which are in qualitative 
agreement with each other. One can simply take the 
temperature at which the difference between l m ax(J,^) 
and, say, l ma x{J = 0, A) is a certain value (or at which 
this difference is a certain percentage of l ma x(J, A)) to be 
T^(I). While the values are arbitrary, it is found that 
the temperature dependencies of l\ (T) and 1% (T) do not 
depend on the chosen values nor on whether one chooses 
the absolute difference or the relative difference. What 
does depend on these factors are the slopes. For this 
reason, it should be emphasized that Fig. [l] is only a ren- 
dition of the phase diagram. The slopes of these lines 
as well as the temperature and current scales will vary 
not only from the less anisotropic materials to the more 
anisotropic material but also from sample to sample. ||] 
The fact that the temperature scale and the current scale 
cannot be determined in our RG analysis contributes to 
this ambiguity. Q 

The temperature dependence of ll{T) has also been 
considered by Jensen and Minnhagen |2| who found it 
to be linear at all currents fcOj unlike the small- cur- 
rent behavior found in Ref. [pL Using mean-field the- 
ory, they found that dl\ / dT oc I/7 where 7 = £ a b / £ c 
is the anisotropy factor. (The relationship between this 
quantity and our parameter for the strength of the in- 
ter layer coupling is 7 oc 1/vA.) This issue could not be 
addressed in our approach for the following reason. When 
the anisotropy is large, « T c and our temperature 
scale, which is linear in t = T/Tkt — 1, is only valid 
for small t. For similar reasons, the dependence of I^iT) 
on 7 could not be determined. For I C (T), no systematic 
dependence on 7 could be identified. 

C. The phase transition at I C {T) 

As mentioned above, a second order phase transition 
of this system occurs at I C (T). It is here that the non- 
analyticity in the correlation length occurs and above 
which y(V) goes to infinity for large I and below which it 
goes to zero. In terms of vortices, we believe that at this 
current, free vortex lines can be spontaneously created 
as opposed to being created by a thermal unbinding of 
vortex loops. This means that the total vorticity of the 
system need not be zero but will fluctuate around zero. 
This leads to the meaning of I^iT) which we see as a 
natural extension of the definition of I C (T). As we have 
stated before, at this current the system starts to be af- 
fected by the current which we have taken to mean that 
it is now large enough for vortex loops to be thermally 
unbound into vortex lines. In other words, there are two 
mechanisms for creating vortex lines, one is by an un- 
binding of vortex loops and the other is by spontaneous 
creation. The latter is possible because the energy of the 
vortex line becomes finite due to screening. In zero cur- 
rent these two mechanisms become possible at the same 



temperature whereas in finite current, they occur at dif- 
ferent temperatures. This is why I C (T) and ll(T) start 
out at the same point in zero current but become separate 
lines at finite current. As mentioned above, Martynovich 
and Artemov fllCfl have studied this system in zero and 
finite current. Without a current, they predict a first- 
order phase transition which turns into a second-order 
phase transition at a finite current. No evidence for this 
has been presented however and in fact there is evidence 
for a second order phase transition at zero-field , pjl ] 

Gr0nbech- Jensen, Dominguez, and Bishop Jllfl have 
arrived at conclusions similar to ours based on Langevin 
simulations of anisotropic 3D, current-driven Josephson 
Junctions arrays. There they calculate two quantities, 
the voltage and the helicity modulus. They find that 
the helicity modulus goes to zero at a temperature above 
which the voltage becomes finite. They conclude that 
the phase transition occurs at a temperature higher than 
that at which vortex loops start to thermally unbind in 
agreement with our results. They also find that I^{T) 
and J C (T) are linear in temperature. 

To summarize this section, Fig. |, where I^(T), I C (T), 
and 1% (T) are plotted, represents the "raw" results of our 
analysis. ll{T) is the value of the current at which the 
system starts to be affected by the current; I C (T) is the 
current at which the phase transition occurs; and Ic(T) is 
the current at which the 3D/2D crossover occurs. These 
results were obtained for a system composed strictly of 
vortex pancakes and neglect the structure and behav- 
ior of the underlying superfluid. Amplitude fluctuations 
of the superfluid become significant near the mean field 
transition temperature T c q and the superfluid is next to 
non-existent making vortices a minor detail. Minnhagen 
p2j has taken this into account and we will see an exam- 
ple of its effect when we consider Eq. (||). (The position 
of T c q(I) relative to T£(I) could not be determined in our 
analysis.) In the next section we will relate our results 
to layered superconductors such as the high- tempera- 
ture superconductors. In particular we will describe the 
expected electric transport properties that are heavily 
influenced by vortices in each region for these materials. 

III. I-V AND R(T) CHARACTERISTICS OF THE 
I-T PHASE DIAGRAM 

In this section, the regions of the I-T phase diagram 
will be discussed and the expected electrical transport 
properties (e.g., TV curves and resistance R(T) mea- 
surements) of each region will be explained. (See Fig. g.) 
Flux transformer measurements and the detection of 
I%(T) are also discussed. 
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FIG. 2. The expected current-voltage characteristics of the 
I-T phase diagram. In Region A, V is expon entially small but 



not zero because, as discussed in Sec. Ill A, vortex loops can 
blow out in a direction perpendicular to the layers. In Region 
B, Eq. ^) holds at small currents (shaded area) but will break 
down at larger currents. Region D has a weak non-ohmic con- 
tribution (Eq. ^)) which enters through the renormalization 
of the transition temperature. 



A. Region A 

This region is defined by where the full correlation 
length and the layered, zero-current correlation length 
are the same. The coincidence of these two quantitites is 
interpreted to mean that in this region the current does 
not have a significant effect on the system. In the original 
paper reporting the I-T phase diagram, the current- 
voltage characteristic was taken to be V(I) = 0. In real- 
ity, it should read V(I) ~ since vortex loops do have a 
exponentially small chance of being unbound in the pres- 
ence of a current Jl],^3| because the current also exerts a 
force on the segments of the loop which lie between and 
parallel to the layers. In a weakly coupled layered system 
this can be neglected and l\ (T) signifies a blowout of the 
loops in a direction predominantly parallel to the planes. 
As one goes to the more isotropic systems, this quantity 
loses its meaning since the loop can blow out in both the 
parallel and perpendicular directions simultaneously. In 
their simulations, Gr0nbech- Jensen et al. do find a 
small number of blow-outs at currents less than I^(T). It 
should be noted that in the two-dimensional limit, this 
region does not exist because I\ (T; 7 = 00) = 0. 



B. Region B 

In this region, the full correlation length and the lay- 
ered, zero-current correlation length differ from one an- 



other signifying that the current has a significant effect on 
the critical behavior. This implies that just above ll{T), 
vortex loops are being thermally unbound in the direc- 
tion parallel to the layers and the I-V relation (Eq. ([j])) 
of Jensen and Minnhagen || should hold. But this equa- 
tion should not be expected to hold in all of Region 
B for a number of reasons. First of all, the equation 
is based upon the formula for rare escapes over a bar- 
rier. p3] This breaks down severely as one increases the 
current and the energy of the most energetic pairs ap- 
proaches that of the barrier height. Secondly, it is based 
upon vortex loop unbinding in the parallel direction (i.e., 
2D unbinding). Close to the transition temperature, 
the system is three-dimensional in its behavior and loop 
blow outs in the perpendicular direction become signif- 
icant. These effects may be represented to first order 
by introducing a current dependence into the parame- 
ters of Eq. (|l|). In other words, to be more strict one 
should write Eq. (|) as V = k{T)I{I-l{T, I)) a( * T D with 
1(T,I ~ Ic( T YL- ^(T)- At lar S er currents, I(T,I) 
goes to zero. |g5[ The behavior of a(T, I) is less clear. 
Our renormalization group analysis showed that a(T, J) 
should be independent of current over a substantial part 
of region B at fixed temperature but this analysis does 
not incorporate either of the effects mentioned above. 
What is certain is that a(T,I) goes to zero at I C (T). 
The approximate region in which Eq. ([[]) is expected to 
hold is represented in Fig. by the shaded region. 



C. Region C 

Region C corresponds to temperatures above T C (I) 
where the full correlation length differs from the 2D 
finite-current correlation length. Because one is above 
the phase transition free vortex lines can be sponta- 
neously created which one typically considers to mean 
that the region is ohmic . However, as we will show 
here and in Section HID, a non-ohmic contribution en- 



ters through the renormalization of the critical temper- 
ature due to the current. There are actually non-ohmic 
contributions from two sources. The first source is of 
course blow-outs of smaller vortex loops which still ex- 
ist above I C (T). The other piece will enter through the 
term describing the spontaneously created free vortex 
lines which is "largely" ohmic. The resistance is propor- 
tional to 1/£(T) 2 but a current dependence enters into 
this through the current-dependent transition tempera- 
ture T C (I). Because our RG analysis does not accurately 
describe the 3D region, the functional dependence of the 
resistance in Region C cannot be determined. The situ- 
ation is better in Region D. 

A curiosity of this region is that its width increases 
with current implying that the 3D region becomes larger 
even though the current tends to decouple the layers. We 
believe that this effect is best looked at by saying that 
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the current has a stronger effect on I C {T) than I^{T). 
It should also be noted that because the full correla- 
tion length differs from the 2D finite-current correlation 
length in this region, the behavior should be of a 3D na- 
ture. At the same time however, the interlayer coupling 
at large length scales is renormalized to zero p6|-p8[] mak- 
ing 2D signatures conceivable. 



D. Region D 

In this region, the full correlation length and the 2D 
finite-current correlation length overlap meaning that 
the behavior of the layered system is 2D-like. One can 
therefore write down the resistance formula for this re- 
gion based on Minnhagen's 2D formula, |^2| R(T) = 
cxp{-[6(r c0 (/) -T)/(T-T K t(/))] 1/2 }, where T c0 is the 
mean-field transition temperature, Tkt is the 2D tran- 
sition temperature and b is a constant. To use this for 
Region D in our the layered system, we must insert the 
current dependences for the two transition temperatures: 

R(T) = exp{-[6(T c0 (/) - T)/(T - T KT {I))]^} , (6) 

This equation is important because it shows how a weak 
current dependence can enter into the resistance formula 
making this region weakly non-ohmic. Eq. (|^) has been 
checked for a YBCO data and it is found to work well 
p9[ for weak currents (< 10 _3 A) using a constant T c0 . It 
should be stressed that Tkt (I) is not the same as T C (T). 
Tkt{I) is the temperature at which it appears that the 
correlation length in the 2D region would diverge, and 
T c (7) is the actual temperature at which it diverges once 



rises rather abruptly from a value zero in contrast to the 
high-temperature side where V s descends to zero but but 
starts to increase before reaching it. 
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the 3D effects take over. As noted above in Section HC 



one should keep in mind that the underlying superfluid 
is also being destroyed in this vicinity and so one must 
also consider normal state effects which can wash out the 
effects of vortices. 



E. Ic(T) and Flux Transformer Geometries 

Here we will show how our results can be used to 
explain measurements made in the flux transformer ge- 
ometry ||(J and explain how one could measure I%(T) 
which is not easily determined from the primarily in- 
plane Eqs. (Q) an d (|6| ). In the flux transformer geom- 
etry arrangement, |30( | a current is injected in the top of 
a sample while a secondary voltage V S (T) is measured 
on the bottom. Typically, the current is held constant 
while the temperature is varied so that horizontal slices 
of the I-T phase diagram can be studied. Such a mea- 
surement has been carried out in zero-field by Wan et 
al. fl3l[ where it was found that near the transition tem- 
perature a peak appeared in V S (T) for various current 
values. On the low-temperature side of the peak, V S (T) 
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FIG. 3. 7 C (T) determined from V B (T) as measured in the 
flux transformer measurements of Ref. f§N- As discussed in 
the text, within the error bars this data is consistent with our 
results that (7) be linear in current 

The behavior of the V s (T) peak can be explained in 
terms of the I-T phase diagram. The onset of the sec- 
ondary voltage on the low-temperature side is T*(J) 
where the loops starts to blow out due to the current. 
Because there is finite interlayer coupling in this region, 
the movement of free vortices in the upper layers is trans- 
lated through to the bottom of the sample. The peak gets 
larger as one goes deeper into Region B because more 
loops are blowing out. However, as one goes into Region 
C, this effect is offset by the gradual layer decoupling and 
so V s (T) begins to decrease and in principle will decrease 
to a value zero at T^(I). If one could neglect amplitude 
fluctuations, V S (T) would decrease to zero at Tj?(I) but 
this was not the case in Ref. |HJ. We have determined 
l\ (T) from the secondary voltage data of Ref. |3lJ using 
an onset voltage of O.OlyitV as our threshold. It is plotted 
in Fig. |§| along with its error bars which were determined 
from the temperature resolution used in the measurement 
of V s (T) . Within the error bars it is possible to say that 
the data is linear in T as found by others using in-plane 
I-V measurements ||-f||| and also predicted in Ref. 
although it is not as convincing as the other data. 

Even though T c 2 (7) could not be determined from the 
data of Ref. |nj this method remains the best hope for 
studying 7,?(T). It may be possible to find a sample in 
which the layers become decoupled before the underlying 
superfluid breaks down and a close inspection of the Wan 
et al. data suggests that (I) might be determined at 
small currents. It should also be noted that the their data 
does manifests the behavior predicted by our results and 
discussed in Section [II C: as the current is increased, 
T c 2 (7) -T^(I) increases. 
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IV. FINITE SIZE EFFECTS 

In this Section, we consider finite size effects which 
make the energy of a single vortex finite and therefore 
make it possible to have free vortices in zero-current 
at temperatures below the critical temperature. The 
presence of these free vortices destroys the Kosterlitz- 
Thouless-Berezinskii (KTB) transition |32],[33| turning it 
into a crossover. This is because the free vortices screen 
the vortex pairs more strongly than do pairs which in 
turn makes vortex pair unbinding possible at lower tem- 
peratures as originally discussed by Kosterlitz and Thou- 
less. J32| Nevertheless, it is still possible to see KTB be- 
havior as we explain below. 

There are two ways in which finite size effects enter. 
In a chargclcss supcrfluid, the bare energy of a 2D, sin- 
gle vortex is In L where L is the size of the system. In 
charged superfluids, the energy of a 2D, single vortex goes 
as In \\jd where Al is the London penetration depth and 
d is the thickness of the film. Usually, however, L and Al 
are large in typical samples and thus do not smear the 
transition in a significant way. Beasley, Mooij, and Or- 
lando p3] along with others ||^] pointed this out in the 
late 1970's. In the presence of a current, the principle 
consequence of a finite L or Al is to contribute an ohmic 
term to the otherwise purely non-linear (below T c ) I-V 
characteristics which can dominate at low currents and 
can wash out the phase transition. 

In the HTSC's there is evidence that Al is small and 
therefore can significantly affect the critical behavior in 
this system. Indeed, linear I-V characteristics at small 
currents have been observed by Paracchini and Romano 
Ji2| , [r3[ on YBCO, by Matsuo et al. on layered conven- 
tional superconductors, jl4| in the author's own anal- 
ysis of I-V measurements from various groups on 
YBCO and BSCCO materials, and by Repaci et al. in 
ultra-thin YBCuO films Jl5|,[l6| who attributed the small 
current ohmic behavior to finite size effects. Repaci et 
al. [ p>[|l6| ] go on to report that low-sensitivity data are 
consistent with a KTB transition. Here, we will address 
finite size effects on the critical behavior and the I-T 
phase diagram. While their effect may ultimately de- 
stroy the phase transition, our position |36| is that KTB 
behavior can still be observed and the structure of the 
I-T phase diagram is reflected. (The effect of a finite L 
on 2D Joscphson Junction arrays was recently considered 
theoretically and in simulations by Simkin and Kosterlitz 

©■) 

As mentioned above, finite size effects introduce an 
ohmic term into the I-V characteristics which are typ- 
ically purely non-linear below T C (I) which can be seen 
through the following simplified derivation. The kinetic 
equation for the formation of free vortices in terms of the 
density of free vortices n F is 

dn F /dt = Tj{T, J) + T FS (T) - k in %, (7) 



where Tj(T, J) is the rate at which free vortices are being 
created by pair unbinding due to thermal activation over 
the barrier made possible by the current, T F s(T) is the 
rate at which free vortices are created due to finite size 
effects, and k\ is a constant. It is a standard derivation 
Q to show that 

Tj{T,J)<x{I~ll{T)) 2 ^ T \ (8) 

which corresponds to the rate at which bound vortex 
pairs are thermally activated over a barrier and as{T) 
corresponds to a bare exponent that must be renormal- 
ized. The rate at which free vortices are spontaneously 
created is proportional to a Boltzman factor: 

T FS {T) cx exp[-E FV /k B T], (9) 

where E F y is the energy of a free vortex. As pointed out 
above, Eq. |§| is valid for rare escapes over a barrier and 
also for when the vortex loop blowouts are in a direction 
parallel to the planes. In that limit, one can take T F s(T) 
to be independent of current. It is Eq. [l] that one uses to 
derive Eq. Ill in the limit of E F y = oo. 
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FIG. 4. The approximate modifications to the I-T phase di- 
agram due to finite-size effects. Note that the original phase 
lines are shown here as dashed lines. Within the solid lines 
is where one would expect the effect of vortex-pair unbinding. 
Outside of this region ohmic behavior is expected although it 
could be exponentially small at lower temperatures and cur- 
rents. Slightly non-ohmic behavior may persist to the right 
on the non-ohmic region as a result of transition temperature 
renormalization. 

Proceeding with Eq. 0, one assumes a steady state so- 
lution and finds 

n F cx y/Tj{T, J)+T FS (T). (10) 

This is readily put in the form of a current-voltage rela- 
tion since R<xn F and V = IR: 
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V = IyJrmiT) + fc(T) 2 [/ - /i . (11) 

At low currents, the first term will dominate and the 
TV curves will be linear but at larger currents it is pos- 
sible for the second term to dominate and for the I-V 
curves to become non-linear. Therefore, by taking into 
account the first term, it is still possible to determine 
J*(T) and a(T). The current I*(T) at which the second 
term dominates introduces a fourth characteristic current 
into the I-T phase diagram which becomes quite different 
in the presence of finite size effects as we discuss below. 
If the temperature dependence of mi(T) and k(T) were 
known, it would be possible to determine the relationship 
between I*(T) and a(T). As it is we have, 

I* (T) = ll(T) + [m 1 (T)/k(T) 2 ] 1 ^ 2a ^ (12) 

which tells us that I* (T) lies above l\ (T) and that their 
difference increases with temperature. (In this last state- 
ment we have made use of some empirical understand- 
ing of m\(T/k(T).) Based on this we can modify the 
I-T phase diagram to include finite- size effects as we 
illustrate in Fig. [|. Within the solid lines is where one 
would expect the non-linear behavior due to vortex pair- 
unbinding. The lower solid line is I*(T) and the upper 
is I C {T). We have not included the effect of the renor- 
malization of I C {T) due to the free vortices here. To a 
first approximation, their effect will be to shift this quan- 
tity to lower temperatures and perhaps to smear it. The 
renormalization of the quantities in Eq. [ll] should also 
be accounted for but here we do not expect any drastic 
changes. Fig. [| is similar to Fig. 8 of Ref. [|l3| but the 
differences are important. Their Iq(T) is defined through 
a phenomenological equation different than Eq. (Q) and 
is therefore distinct from (T) and I*(T). Furthermore, 
the ohmic region and non-ohmic regions have different 
shapes in the two figures. 

V. EFFECTS OF DISORDER 



The essence of the replica technique is to average the 
disorder over n replicas of the system which is done via 
the identity [InZ] = lim T ,_ >0 / dV d P(V d ){Z n -l)/n where 
Z is the partition function for the layered vortex gas and 
[...] represents an overage over the probability distribu- 
tion of the disorder. P{V d ) is a Gaussian distribution 
with zero mean value and short ranged spatial correla- 
tions: [Vd(r)Vd(r')] ~ A8(\r — r'|) where A represents 
the strength of the disorder. The effect of this averaging- 
is to introduce an interaction between vortices in differ- 
ent replicas which includes the standard logarithmic in- 
teraction in addition to a short ranged disorder induced 
interaction which is proportional to — /3AS(\r — r'|) as de- 
scribed in more detail in Ref. [[l3| After this averaging, 
the model is similar to that of Fischer except the ad- 
ditional disorder induced interaction is long ranged there. 

We have done a renormalization group analysis on 
the new partition function which consists of two steps: 
a course graining and a rescaling. While the integra- 
tions for the course graining step cannot be done exactly, 
one can show that the resulting terms have no logarith- 
mic or linear dependencies and so will not renormalize 
the interaction strength parameters p 2 and A. To do 
this integral, one approximates the delta function by the 
short-ranged expression, jl3| exp[— r 2 /£ 2 ]. The rescaling 
steps alone result in the following recursion relation for 
the disorder strength: dA/de = — A which is immedi- 
ately seen by considering —/3A6(R) with the substitu- 
tion R' = R/(l + dr/r). Since it is unlikely that coarse- 
graining effects will contribute to making the disorder 
more important, this parameter is irrelevant. 

Because the disorder does not affect the recursion re- 
lations of the system parameters p 2 , A and y and the 
disorder parameter A is irrelevant, our results show that 
a weak, random, one-body potential has no effect on the 
critical behavior of the system in agreement with the re- 
sults of Ref. |4^] . This is expected since even with one 
vortex in a pair pinned, the pair can still unbind because 
of the "freedom" of the other vortex. 



In this section, the effect of a quenched, random one- 
body potential Vd(R) on the system will be considered 
using the replica technique. |pq| The primary effect of 
this type of disorder is to pin a vortex, but various types 
of disorder on 2D systems have been considered in the 
past. Rubinstein et al. J3?J and Korshunov j!o| have 
studied the effects of a random configuration of quenched 
vortex pairs. Jose has investigated the 2D ferromag- 
netic planar model where the coupling is taken to be a 
random variable while Fischer jl2| has studied a system 
where the disorder is introduced via the fugacity. The ef- 
fect of a quenched, random one-body potential on a one- 
component vortex gas (i.e., no antivortices) has been con- 
sidered by Menon and Dasgupta Q] and we will follow 
their notation here. We will present the calculation for 
the layered case which is also applicable to the 2D case. 



VI. DISCUSSION AND SUMMARY 

Numerous aspects of the current-temperature phase di- 
agram have been considered here including its derivation, 
the current- voltage characteristics of each region, proper- 
ties of the phase lines, and the use of the flux transformer 
geometry to study I 2 (T). We have also considered finite 
size effects and derived a current-voltage equation for this 
case. Finally, the effect of a weak, random, one body po- 
tential was found to be weak. An interesting consequence 
of the current dependence of the critical temperature is to 
make the system slightly non-ohmic above T c . There re- 
main many interesting aspects of the I-T phase diagram 
to consider including the dependencies of the character- 
istic currents I^iT), I C (T), and I 2 (T) on the anisotropy 
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factor 7. 
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